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ABSTRACT 

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for revisions 
in the picture of particle acceleration. We make use of the recently established scaling of slow and 
fast MHD modes in strong and weak MHD turbulence to provide a systematic study of particle 
acceleration in magnetic pressure (low-/3) and gaseous pressure (high-/?) dominated plasmas. We 
consider the acceleration by large scale compressions in both slow and fast particle diffusion limits. 
We compare the results with the acceleration rate that arises from resonance scattering and Transit- 
Time Damping (TTD). We establish that fast modes accelerate particles more efficiently than slow 
modes. We find that particle acceleration by pitch-angle scattering and TTD dominates acceleration 
by slow or fast modes when the spatial diffusion rate is small. When the rate of spatial diffusion of 
particles is high, we establish an enhancement of the efficiency of particle acceleration by slow and 
fast modes in weak turbulence. We show that highly supersonic turbulence is an efficient agent for 
particle acceleration. We find that even incompressible turbulence can accelerate particles on the 
scales comparable with the particle mean free path. 
Subject headings: MHD — turbulence — acceleration of particles 



1. INTRODUCTION 

MHD turbulence is an important agent for particle 
acceleration as was pointed first by Fermi (1949) and 
later was discussed by many other authors (see Chan- 
dran & Maron 2004b and references therein). Second 
order Fermi acceleration by MHD turbulence was ap- 
pealed for acceleration of particles in many astrophysical 
environments, e.g. Solar wind, Solar flares, the intraclus- 
ter medium, gamma-ray bursts (see Schlickeiser & Miller 
1998; Chandran 2003; Petrosian & Liu 2004). Naturally, 
properties of MHD turbulence (see Shebalin, Matthaeus, 
& Montgomery 1983; Higdon 1984; Montgomery, Brown, 
& Matthaeus 1987; Shebalin & Montgomery 1988; Zank 
& Matthaeus 1992; Cho & Lazarian 2005 and references 
therein) are essential for understanding the acceleration 
mechanisms. 

For most part of the paper we use the model of strong 
MHD turbulence initiated by a pioneering work by Gol- 
dreich & Sridhar (1995, henceforth GS95). GS95 dealt 
with incompressible MHD turbulence and showed Alfven 
and pseudo- Alfven modes follow the scale-dependent 

anisotropy of Z|| ~ L 1 / 3 / 2 / 3 , where Z|| is the size of the 
eddy along the local mean magnetic field, lj_ that of the 
eddy perpendicular to it, and L the outer scale of turbu- 
lence. Lithwick & Goldreich (2001) conjectured that this 
scaling of incompressible modes is also true for Alfven 
modes and slow modes in the presence of compressibil- 
ity. In Cho & Lazarian (2002; 2003, henceforth CL03) 
we provided arguments justifying the rational for consid- 
ering separately the evolution of slow, fast and Alfven 
mode cascades. Our numerical simulations 1 verified that 
Alfven and slow mode velocity fluctuations indeed show 

1 Higher resolution numerics in Cho & Lazarian (2004) confirms 
our earlier results. 



GS95 scaling, while fast modes exhibit isotropy in both 
gas-pressure (high (3) and magnetic-pressure (low j3) 
dominated plasmas (see review by Cho & Lazarian 2005). 

MHD turbulence may also be in weak regime (Galtier 
et al. 2000). In weak MHD turbulence the cascade in- 
creases only the perpendicular wavenumber and does not 
change the parallel one. The inertial range of the weak 
turbulence in most astrophysical circumstances is lim- 
ited. As the turbulence gets more and more anisotropic 
it transfers into the strong MHD turbulence 2 . In spite of 
the limited inertial range, the weak turbulence may be 
important for cosmic ray acceleration. 

The new paradigm of MHD turbulence has substan- 
tially changed our understanding of energetic particle- 
turbulence interactions via gyroresonance and the 
Transit-Time Damping (TTD) (Chandran 2000; Yan & 
Lazarian 2002, 2004; Farmer & Goldreich 2004). How- 
ever, these two processes do not exhaust all the relevant 
interactions. For instance, the acceleration of cosmic rays 
by the large scale compressible motions was described in 
the literature rather long ago (see Ptuskin 1988, hence- 
forth P88). Attempts of applying the GS95 scaling to 
the stochastic large-scale MHD acceleration of cosmic 
rays were described in some recent publications (Chan- 

2 The names strong and weak may be misleading therefore. 
Strong MHD turbulence may be of very small amplitude. Its char- 
acteristic property is a so-called critical balance, which is the bal- 
ance between the time scales of perpendicular and parallel motions, 
i.e. l±/vj_ ~ where va = Bo / y/4np is Alfven speed. Here 

Bq is the strength of the mean magnetic field and p is density. The 
weak turbulence decreases the perpendicular scale of the eddies 
while keeping the parallel scale Zm intact. This process can proceed 
till the critical balance condition is reached. After that the turbu- 
lence is in the strong MHD regime and both l± and lu change in 
a way to preserve the critical balance. These issues are discussed, 
for instance, in Cho, Lazarian & Vishniac (2003). 
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dran 2003, henceforth C03; Chandran & Maron 2004ab). 
The results were applied to particle acceleration in Solar 
Flares and clusters of galaxies. 

These studies raise a number of questions. For in- 
stance, it was found in Yan & Lazarian (2002, hence- 
forth YL02) that fast modes dominate scattering of cos- 
mic rays if turbulent energy is injected at a large scale. 
Do the slow modes, nevertheless, dominate the particle 
acceleration in Solar flares? What is the relative contri- 
bution of large scale compressions and small scale scatter- 
ing to the particle acceleration? Do the results change if 
the turbulence is weak (see Galtier et al. 2000)? Is com- 
pressibility a necessary requirement for acceleration? Is 
there a possibility to reconcile results by P88 and C03? 
These and other questions induce us to undertake a de- 
tailed study of particle acceleration by MHD turbulence. 

We consider particle acceleration by both strong and 
weak MHD turbulence. We consider fast and slow modes, 
highly supersonic motions, effects of betatron electric 
field. We compare the rates of acceleration by differ- 
ent non-resonant mechanisms and the acceleration rates 
arising from pitch- angle scattering and TTD. 

For the most part of the paper, we consider a fluid 
threaded by a strong mean magnetic field Bo- But the 
results may easily be generalized if local magnetic field 
is used instead of B (see CL03). 

We start with introducing of the basics of the non- 
resonant acceleration by compressible turbulence in §2. 
In the same section we discuss particle acceleration in 
acoustic turbulence. We consider particle acceleration 
by slow modes (§3) and fast modes (§4), highly super- 
sonic turbulence (§5), incompressible turbulence (§6), 
and weak turbulence (§7). We discuss gyroresonance and 
TTD acceleration in §8 and compare the efficiency of dif- 
ferent mechanisms in §9. We provide discussion in §10 
and summary in §11. 



2. ACCELERATION BY COMPRESSIBLE MOTIONS: 
HYDRODYNAMIC APPROACH 

2.1. Basics of non-resonant acceleration 

The acceleration by large-scale motions depends on 
compression created by different modes. The traditional 
method of estimating D p is based on calculation of the 
degree of fluid compression within a hydrodynamic ap- 
proach. It is implicitly assumed that the motions in ques- 
tion are much larger than a particle mean free path, the 
magnetic field moves together with the fluid and particles 
are moved together with the magnetic field. 

When a particle moves inside an eddy of size I, the 
change of particle's momentum over a time At is 



The change of momentum depends on compression 



3. 



Ap ~ (dp/dt)At, 



(1) 



where p is the momentum of the particle. The momen- 
tum diffusion coefficient is 



dp 



777 ~ ~P~ 

dt Co 



(3) 



where a constant Co depends on the geometry of compres- 
sion. For example, cq = 3 when compression is isotropic 
(see P88). Here, for the sake of simplicity, we set cq = 1. 
Substituting eq. J2J into eq. we obtain 



D p ~p 2 {V -v t ) 2 At. 



(4) 



Therefore we can determine D p when we know V • v; and 
At. 

2.2. Wave-like turbulence 

An important work for stochastic acceleration is P88, 
that implicitly considers acoustic turbulence 4 . Ptuskin 
compared two time scales on the outer scale of turbu- 
lence - wave period (t wave ~ L/a) and diffusion time 
(tdiff ~ L 2 /D), where a is the sound speed, D the spa- 
tial diffusion coefficient, and L the outer scale of turbu- 
lence. 

Fast diffusion limit 
When diffusion time, tdiff, is smaller than wave period, 
twave, we can use At = tdiff- P88 considered isotropic 
turbulence arising from interacting sound waves. For the 
outer scale eddies, since V • vj, ~ Vl/L, we get 

D p ~p 2 (v 2 /L 2 )(L 2 /D)~p 2 v 2 L /D, (5) 

where vl is the rms velocity (P88). 
It is straightforward to calculate D p for smaller eddies: 



D, 



p 2 v A Lv A ( vi 



D 



va 



(6) 



for eddies on scale I < L, which is smaller than D p for 
outer scale eddies. When we want to calculate con- 
tributions from all scale eddies, we need to calculate 

jD p (di/iy. 



D p {dl/l) ~ D p>L I (l/L) 2m {dl/l) ~ D p . L /2m, (7) 



where D p l is given in eq. 10 and we assumed vi ~ 
v L {l/L) m '. Zakharov & Sagdeev (1970) claimed that 
m = 1/4, or equivalently E(k) cx fc -3 / 2 for weak acoustic 
turbulence. The result is very close to the one in eq. 

Slow diffusion limit 
On the other hand, for tdiff > twave, P88 considered 
outer scale eddies and obtained 



D p ~p 2 (Dv 2 L /a 2 L 2 ) 

„2„,2 / 2\ 



~ (P v L/a )/t L ,diff 
(p 2 vl/D)(t L , wave /t L: diff) 2 



(8) 
(9) 
(10) 



D p ~ (Ap) 2 /At ~ {dp/dt) 2 At. 



(2) 



3 We will discuss why the change of momentum depends on V • v 
in 86.31 We can understand it as follows. First, consider compres- 
sion only in the directions perpendicular to the mean magnetic 
field. This will result in compression of magnetic field lines, which 
increases the perpendicular momentum of the particle. Second, 
consider compression only in the directions parallel to the mean 
magnetic field. This will cause adiabatic heating and increase the 
parallel momentum of the particle. Overall, the rate of momentum 
change depends on the rate of compression. 

4 Formally this corresponds to weak isotropic turbulence caused 
by interacting sound waves. In this paper (see §7), weak turbulence 
means that the wave- wave interaction is so weak that wave packets 
can travel distances much larger than their typical size. 
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which goes to zero when D — > 0. This relation can 
be understood as follows. Consider a particle inside a 
region of size L that is compressed/expanded by large 
scale wave packets in weak turbulence. When the parti- 
cle diffuses out of the region, the change of momentum 
is Ap ~ (dp/dt)t L ^ wave ~ (pv L /L)(L/a) ~ pv L ja. How- 
ever, since a particle stays inside the region for approxi- 
mately tL t diff, we take At ~ tL,diff- Substituting these 
Ap and At into D p ~ (Ap) 2 / At, we obtain the result in 
eq. ©. 

According to P88 original formula in eq. 1(8} , particle 
acceleration is inefficient in slow diffusion limit. How- 
ever, note that Ptuskin's formula in eq. JHJ is derived 
for large scale eddies (or, to be precise, large scale com- 
pressive motions). It is easy to show that small scale 
eddies can provide efficient acceleration. Let us consider 
eq. (|lt)p . For large scale eddies, the ratio tz^ave/tL^diff 
is very small in slow diffusion limit. As we move down 
to smaller scale eddies, ti >wave scales as l/a and as 
l 2 /D. The ratio ti_ wave /ti^if / becomes larger as the scale 
decreases. Therefore, as we move down to smaller scales, 
we will have higher acceleration efficiencies. When the 
inertial range is wide enough, we will ultimately reach 
the scale l c on which t^ waV e ~ ti,diff- On the scale Z c , 
the momentum diffusion coefficient is 



D p ~ p 2 v?/D, 
where vi is the velocity at scale l c 



(11) 



3. ACCELERATION BY COMPRESSIBLE MOTIONS: SLOW 

MODES 

Here we consider acceleration by slow modes in mod- 
erately compressible turbulence for low Mach numbers. 
As above, we discuss both slow and fast diffusion lim- 
its. We assume turbulence is strong on the outer scale 
(i.e. vl ~ va)- In this case, the outer scale eddies are 
isotropic. However, smaller eddies are elongated along 

the local mean magnetic field and we have h ~ L 1 / 3 ? 2 / 3 , 
which states that anisotropic eddy shape is more pro- 
nounced on smaller scales. We consider separately gas 
pressure (high-/?) and magnetic pressure (low-/?) domi- 
nated plasmas. As in C03, we assume that spatial diffu- 
sion of particle occurs due to pitch-angle scattering. 

3.1. Fast diffusion limit in low-(3 

Eq. (@J is a universal formula for D p and we can use it 
for slow modes. Evaluation of eq. (@J requires V • vi )S iow 
and At. Since we consider the case of tL,diff < th,eddy, 
where tL,eddy ~ L/va is the eddy turnover time on the 
outer scale, particles diffuse out of an eddy before it is 
randomized. Therefore, we take the diffusion time for 
At: At ~ L 2 /L>||, where we used Dm because particle 
diffusion occurs along the magnetic field lines. For slow 
modes in low-/? limit, we have V • vl.sIow ~ Vl,sIow/L\\ 
(see Appendix B). When turbulence is isotropic on the 
outer scale, as assumed in this paper, we have Lu ~ L. 
The momentum diffusion coefficient is 



r~iSlow 
P 



p 2 (v L ,sio W /L) 2 (L 2 /D ll ) 

p 2 V A Lv A I 



L Dn 



VA 



(12) 



which is very close to eq. (0 for acoustic turbulence. It 
is also straightforward to calculate D p for smaller eddies: 



D 



p,i 



p 2 v A Lv A ( visi 



L D 



va 



(13) 



for eddies on scale I < L, which is smaller than D p for 
outer scale eddies. All these results are similar to those 
for acoustic turbulence. 

The only difference stems from the fact that when par- 
ticles diffuse along magnetic field lines, they can reenter 
the same outer scale eddy (Chandran & Maron 2004b). 
This results in a somewhat uncertain factor N which 
ranges from unity to several (see discussion in the Ap- 
pendix A). This results in 



D„ 



p 2 V A Lv A ( Vi.si 



L Dm 



VA 



N 



(14) 



for eddies on scale I. 



3.2. Slow diffusion limit in low-ft 

Consider d min VA <C D\\ <C Lv A , where d mm is ei- 
ther the parallel size of eddies at the dissipation scale 
or the mean free path of the cosmic rays, whichever is 
larger. The particle diffusion time is smaller than the 
eddy turnover time at the outer scale L and larger than 
that on the scale ~ d m i n . When Dm -C Lva, particle dif- 
fusion time th,d%ff is larger than the eddy turnover time 
tL,eddy at the outer scale L. In this case, for eddies whose 
parallel size is larger than D^/va, particles are confined 
within the eddies until the eddies are randomized. 

For this regime C03 obtained (see Appendix A) 



D, 



P 2 va ( v lmo w x 2 

VA 



for a scale I and 



D, 



p 2 VA 



VL,slow 
VA 



hi 



\ D \\ 



(15) 



(16) 



for action of all eddies. We observe, that each eddy whose 
parallel size is between Du/va and L makes an equal 
contribution to D p . More importantly, the efficiency of 
acceleration becomes large as Dn — > 0. 

On the contrary, when we apply the approach in i J2.2l 
(see discussion below eq. {HIJO to slow modes, we ob- 



tain' 3 D 



V 



v I,sIow/ v aH\)- As in acoustic turbu- 
lence, we have D p — > as D\\ — > 0. Note that the appar- 
ent disagreement between P88 and C03 does not stem 
from isotropy/anisotropy of turbulence. 

C03 assumes that the random walk argument is appli- 
cable to the momentum diffusion. This is a somewhat 
subtle issue. For instance, this is not true when diffusion 
time is much longer than the eddy turnover time. Let us 
consider an extreme case of Dm — > 0. In this case, the 

5 When the particle diffuses out of the eddy, the change of mo- 
mentum is Ap ~ {dp/dt)t l>eddy ~ (pvi tS i ow /l\\)(l\\/v A )- The time 
that the particle spends inside the eddy is ~ t\ MS Si so that we take 
At ~ ti,diff- Substituting these Ap and At into D p ~ (Ap) 2 / At, 



we obtain D p ~ (p 2 vl sloul /v A )(D [{ /l 2 ) - 
which is, in many aspects, similar to eq. JSJ 
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particle and the fluid element move together. The den- 
sity of the fluid element follows the continuity equation: 



dp _ 
_ = „ v . V p 



Rearranging this, we get 



f Dp 
~p~Dt 



Dlnp 
Dt 



-V-v, 



(17) 



(18) 



where D/Dt = d/dt + v • V. We use D/Dt because 
particle is moving with the fluid in the limit of Du — > 0. 
Since 

I Dp D\np , N 

_V-v = y - = £ 19 

p Dt Dt ' v 7 

where p is the momentum of the particle, hip and In p 
behave similarly. Since hip cannot go to infinity, hip 
cannot increase indefinitely. See Appendix B for more 
discussion about the density fluctuations in compressible 
fluids. Therefore, random walk may not be a good ap- 
proximation for slow diffusion limit. This is a kind of 
suppression of random walk in extremely slow diffusion 
limit. 

Suggested Solution — The apparent disagreement of 
earlier works may be reconciled when we consider tur- 
bulent diffusion. Here we mostly talk about diffusion in 
the directions perpendicular to the mean field. In weak 
turbulence turbulent diffusion is negligible and mobil- 
ity of particles relative to the rest of the fluid relies on 
Du only. Therefore we can apply our argument in the 
previous subsection and D p goes to zero as Dm goes to 
zero. However, in strong MHD turbulence, turbulent dif- 
fusion can be efficient and a particle can move to another 
eddy within one eddy turnover time. Recent numerical 
simulations (Cho et al. 2003) support the idea that tur- 
bulent diffusion is indeed efficient. However, diffusion 
in MHD turbulence is a complicated issue that critically 
depends on mobility of magnetic field lines and, hence, 
the efficiency of magnetic reconnection in turbulent en- 
vironments (see Lazarian & Vishniac 1999). Therefore, 
we will not discuss this further here. 

If turbulent diffusion is fast, particles can move to an- 
other independent eddy within one eddy turnover time 
and random walk-like behavior will be fully recovered. 
The result will be the same as the one in eq. (I16|) . 

If turbulent diffusion is slow we can still have an effi- 
cient acceleration (see a similar argument for hydrody- 
namic wave-like turbulence in M2.2|) . Indeed, for slow 
turbulent diffusion, the eddies on the energy injection 
scale induce smaller D p than the estimate in eq. (|15(l . 
We expect that the suppression of random walk will be 
less severe for smaller scale eddies because the inequal- 
ity, D\\ -C 1\\va becomes milder for smaller eddies. We 
expect that when 

Du ~ hv A (20) 



on a (parallel) scale 

h,\\ ~ D\\/va ~ ImfpVptl/VA, 



(21) 



random walk behavior is fully recovered. Here l m f p is 
the mean free path by pitch-angle scattering. Therefore, 
the momentum diffusion coefficient is largest for eddies 
with the parallel size ~ l c u. Note that tdiff ~ t e ddy on 
the scale i c ii, where tdiff is the diffusion time and t e dd y 



is the eddy turnover time. On this scale, the change in 
momentum during one eddy turnover time is 



y dt 



P 



VLslo 



Vl.slo 



In / \ r \ I ' va 

where we used dp/dt ~ pV -vi^iow ~ pvi lS iow/l\\ and At 
1\\/va (see Appendix A). The net momentum diffusi 
coefficient is similar to eq. (I15f) : 

D S J° W ~ (p 2 vlloJvAl\\) 
^L,slow 



(22) 



p 2 v A 



(23) 
(24) 



L \ v A 

where we used vi lS i ow ~ v L ^iow(l/L) 1/3 , and l\\ ~ 
L i/3;2/3 (GS95). Note that this result becomes identical 
to that in equation (|10f> when we use the fact Du ~ 1\\Va 
(or, U >e ddy ~ ti,diff)- One can still characterize the 
particle acceleration arising from slow modes at large 
scales using eq. (|10ll . but this contribution will be sub- 
dominant to that by smaller scales given by eq. (|24ll . 

To summarize, the momentum diffusion coefficient in 
slow diffusion limit depends on turbulent diffusion. We 
can write 

2 / \ 2 

P V A I VL,slo 



d: 



L 



VA 



Q 



TD, 



(25) 



where Qtd ~ 1 when turbulent diffusion is slow (see eq. 
Q24[0 and Qtd ~ ^-(Lva/Du) when it is fast (see eq. 

03). 

3.3. Slow and fast diffusion in high-(3 limit 

So far in this section, we considered acceleration by 
slow modes in low-/? plasmas: (3 < 0(1). In this case, 
V • Vi !S i ow ~ Vi t gi ow /ln. However, when j3 — > oo, we have 

V ' Vl slow ~ Vl s i ow 

W\\) ( 26 ) 
(see Appendix A). We can use the same At: At ~ 1\\/va- 
Therefore, in high-/? limit, eq. (|12|l becomes 



D 



slow 



P 2 va Lva ( Vl 



L D \v Ai 
for fast diffusion limit and eq. (|24|l becomes 



j~\slow 
V 



p 2 v A 



^L.slow 



{TD 



(27) 



(28) 



L \ va 

for slow diffusion limit. See eq. (|25|l for definition of 
Qtd- Evidently due to the (3~ 2 factor the acceleration by 
slow modes in high beta plasma is suppressed compared 
to the case of low beta plasma. 

4. ACCELERATION BY COMPRESSIBLE MOTIONS: FAST 

MODES 

It is easy to see that arguments in the previous sec- 
tion (SjjJJ are applicable to the acceleration of cosmic rays 
by fast modes, which are shown to be isotropic (Cho 
& Lazarian 2002). As we will see, fast modes, which 
are essentially compressions of magnetic field, can be 
very efficient. Especially in low-/? plasmas, compression 
occurs in the directions perpendicular to the magnetic 
field and there is no adiabatic loss/gain in parallel direc- 
tions. Therefore the acceleration mechanism in H6. II will 
be most relevant. Here we consider acceleration by fast 
modes in moderately compressible turbulence, in which 
shock formation is marginal. 
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4.1. Fast diffusion limit 

In fast diffusion limit the acceleration by fast and slow 
modes is very similar. Indeed, consider Dn 3> Lcf, where 
Cf is the propagation speed of the fast wave and depends 
on the plasma /?: cj ~ va for low-/? and Cf ~ a (=sound 
speed) for high-/?. In this limit, particle diffusion is so 
fast that particles diffuse out of even largest fast mode 
wave packets, or eddies, before the eddies complete one 
oscillation. When a particle diffuse through an isotropic 
fast-mode eddy of size I, the change in momentum is 

I 2 Vl.fn.rtl 



Ap~^At 
dt 



PVlJa st 
I 



P 



Vl, fas t 1 - 



(29) 



where we used V • Vij as t ~ v ijast H and At ~ l 2 /D t 
Therefore, 

D p , l ~{Ap) 2 /M~ P 2 v 2 lJa jD h (30) 

which is most efficient on the outer scale L. 6 
Setting I = L, we obtain 

D f ast _ P 2 VA Lc f V A { VLJast V ^ 
P L D l{ Cf \ V A J ' 

If VLj as t ~ va on the outer scale of turbulence and when 
Dn > Lcf, we have 



D fast _ P 2v A ( LCf_^ < P 2 V A 



(32) 



L \ Dn J 
for low-/? plasmas and 

„ tvA ^ r 1/2 < 1/2 (33) 

for high-/? plasmas. Note that the result in eq. I|31|l 
is very similar to that for Ptuskin's acoustic turbulence 
(eq. (©). This is not so surprising because fast modes in 
MHD are similar to weakly interacting isotropic acoustic 
waves (see Cho & Lazarian 2002). 

Similarly to the case of slow modes, factor N that char- 
acterizes the probability of the return of the particles 
back to the accelerating eddy is applicable. Therefore 
eqs. 1)3 1[L 132(1 and i|33|) should be multiplied by this 
factor. With this addition eq. lpTT|) gets analogous to 
the eq. (fT4"jl . Although /? factors appear both for slow 
and fast modes in eqs. (123 and respectively, they 
have different origins. For slow modes this corresponds 
to the suppression of efficiency, while for fast modes the 
/?-factor reflects the fact that the injection velocity is va- 

6 This statement remains true even the case particles re-enter 
the same eddy multiple times. As in Chandran & Maron (2004b), 
the re-entry factor is given by N ~ min{^J D\\T Lrand /l, M lifaat } , 
where t\ ra nd is the eddy randomization time of the fast modes 
on scale I. Although r; ran( j is uncertain, it is certain that 
ti rand > tl e< idy because fast modes cascade may be slower than 
Alfvenic cascade, where ti^ddy is the eddy turnover time of Alfvenic 
turbulence. In this case, the first term inside the parenthesis scales 
with l q , where q < l~ 2 / 3 . The factor Mij ast is roughly ~ z 3 /l, 
where z s is the distance along magnetic field lines over which two 
adjacent magnetic field lines get separated by the distance I. In the 
presence of strong Alfvenic turbulence, z s ~ L 1 / 3 / 2 / 3 . Therefore, 
the second term inside the parenthesis scales with l~ 1//3 . There- 
fore, the second term is smaller and TV oc Z -1 / 3 . The product of 
Dp in eq. I3UI and N scales with f 2 ^ ast '~ 1//3 - If fast modes follow 



scaling of acoustic turbulence, v? , , 

& ' l. fast 



: Z 1 / 2 . If fast modes follow 

l , J <1SI 

Kolmogorov scaling, u- 

decreases as I decreases. Therefore, the outer scale renders largest 
contribution to D 



, , oc Z 2 / 3 . For either case, v? f .1 x / 3 

L,jast ' t,jast 



4.2. Slow diffusion limit 

In this subsection, we consider d min Cf <ti Du -C Lcf, 
where d m i n is either the size of eddies at the dissipation 
scale or the mean free path of the cosmic rays, whichever 
is larger. Particle diffusion time is smaller than the wave 
period at the outer scale L and larger than that on the 

SCale ^ ^mm • 

In principle, the diffusion coefficient for motions at a 
scale I > l c , where 



lmf P V p tl/c f ~ | £[0-1/2 h^/J 



low-/? 



(34) 



(see eq. (|2*TJ) for l c j), can be obtained using the approach 
in assuming that the coherence time for a diffusing 
particle is equal to the period of the fast wave, rather 
than the cascading time. 

When Dn <C Lcf, we have l c ~ Dn/cf <C L. However, 
we do not need to calculate D p for all scales between l c 
and L. As we will see later in eq. (|36f) . the value of D p 
is largest on the smallest scale l c . Therefore it suffices to 
calculate D p on the scale of ~ l c . 

The momentum change Ap during the diffusion time, 
which is the same as the wave period on the scale of Z c , 
is 

Ap - p(vijast/lc)(lc/cf). (35) 

The momentum diffusion coefficient, D p , is 

fast 



D 



p,l 



{ApY/At 



[p^-fliljcf) 



(P 2v tfast/ C f l c 



- (p 2 v 2 LJast /c f l c )(l c /L) 2m 
(jp 2 v 2 LJa Jc f L)(l c /L) 2 ™- 1 



P 2 VA V L,fast V A 



L 



(lc/L) 



2m -1 



(36) 

(37) 
(38) 

(39) 



where we assumed vij as t ~ VLj ast {l/L) m . Note that 
m = 1/3 for Kolmogorov turbulence and m = 1/4 for 
acoustic turbulence (see Cho & Lazarian 2002). 

The difference between the slow diffusion limit for fast 
and slow modes is evident from the comparison of eq. 
(|39|l and eq. I|23|l . The combination that enters eq. I|39|l 
v t fast 1 1 increases as I decreases, while the combination 
that enters eq. v 2 slow /l\\ is independent of I. Conse- 
quently, the rate of turbulent diffusion does not change 
much for fast modes. Indeed, when turbulent diffusion 
is fast, we do not have suppression of random walk for 
large scale eddies. Therefore, we may need to consider 
all scale eddies between L and l c . However, this does not 
really matter, if most contribution comes from the scales 
near the scale l c for which the suppression is marginal 
even for slow diffusion limit. 

For low-/? plasmas, ct ~ va and fast modes give 



Jjfo-s 



p 2 v A 



P 2 va 
L 

VLjast 
VA 



VLjast 
VA 



L 



1- 2m 



Lva 



V 



-2m 



(40) 



which is larger than that of slow modes if l c <C L. Note 
that l — 2m > for m < 1/2. For high-/? plasmas, ct ~ a 



(> 
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(=sound speed) and fast modes give 

■1 



l-2m 



P 2 VA ( V LJast \ 2 



Lv A 



1— 2m 



, (41) 



L \ Va J \lmfpVptlJ 

which is larger than that of slow modes if l c <C L. 

5. ACCELERATION IN COMPRESSIBLE HIGHLY 
SUPERSONIC TURBULENCE 

When turbulence is highly supersonic, the treatment 
of V • v should be different. In subsonic case, when we 
consider motion on scale I, the general treatment is that 
Ap ~ (dp /dt) At ~ (V-v)At ~ Atvi/l ~ constant, where 
we used At ~ i/wj. However, in the supersonic case, 
shocks compresses gas and form high density regions (see, 
for example, Padoan, Nordlund, & Jones 1997; Padoan 
& Nordlund 1999; Beresnyak, Lazarian & Cho 2005; Kim 
& Ryu 2005). 

The maximum compression for isothermal hydrody- 
namic turbulence scales as ~ Mg, where M s is the sonic 
Mach number (Padoan et al. 1997). We can derive this 
by equating turbulence pressure, ~ pv\, and gas pressure 
in the high density regions, PmaxC 2 , where p is the av- 
erage density and p ma x is the density of the compressed 
regions. Here we assume isothermal gas. The result is 

Pm a X /p~vl/c 2 s =Ml (42) 

This scaling is valid for hydrodynamic turbulence. 

It is not clear whether or not such a simple scaling 
exists for MHD turbulence (see discussions in Ostriker, 
Stone, k, Gammie 2001; Padoan & Nordlund 2002). In 
MHD turbulence, magnetic pressure provides resistance 
for compression, which results in less compression. Since 
slow modes are almost parallel to the local mean field 
directions, magnetic pressure may not provide a signif- 
icant pressure when slow modes form shocks. There- 
fore, we expect that the maximum compression in MHD 
turbulence depends on the sonic Mach number even for 
strongly magnetized case. 

Our numerical simulations support this idea. We use 
data cubes obtained from direct compressible MHD tur- 
bulence simulations with 216 3 cells. These simulation are 
described in CL03. We used 3 data cubes with different 
sonic Mach numbers: M s ~ 0.5, 2.2, and 7. In all sim- 
ulations, we fix the r.m.s. velocity and the strength of 
the mean magnetic field. The Alfven speed of the mean 
magnetic field is very close to the r.m.s. velocity of tur- 
bulence. To see how much compression is achieved, we 
plot the relation between volume filling factor and frac- 
tion of mass contained inside the volume for high density 
regions (see left panel of Fig. ^ . The figure shows that 
higher sonic Mach number fluids are compressed more, 
which is not so trivial because the fluid is strongly mag- 
netized. Let use focus on diamond symbols in the figure. 
For example, the diamond on the dashed curve in the left 
panel says that 25% of mass occupies 2.7% of the total 
volume. 

It is not clear how to derive exact scaling relations from 
Figure n However, we can derive at least two facts from 
the figure. First, as we mentioned earlier, the compres- 
sion rate scales with even for strongly magnetized 
medium. Second, the right panel of Fig. show that 



the compression rates in MHD and hydrodynamic tur- 
bulence scale very similarly. The dashed line in the right 
panel is for hydrodynamic turbulence. Comparing the 
solid line and the dashed line on the right panel, we note 
that the slopes are very similar. It is natural to assume 
that the compression rate in hydrodynamic turbulence 
scales with M 2 . Therefore, we can say that the compres- 
sion rate, defined as the volume filling after compression 
divided by that before compression, in MHD scales with 
M^, where 7 is slightly smaller than 2. 

Here we assume that the maximum compression in a 
strongly magnetized MHD turbulence scales with M] . 
We also assume that spatial diffusion of particle is slow. 
Then we have Ap/p ~ pmax I P ~ The momentum 

diffusion coefficient is 



y\super sonic 
V 



(Ap) 2 /At~p 2 M^(v L /L), 



P VA VL 

L v A 



M 2 S \ (43) 



where we consider only the largest scales because smaller 
scales are less efficient. If vl ~ va, this diffusion co- 
efficient is ~ M s 27 times larger than the one in C03. 
Note that this is a second order Fermi process associ- 
ated with shocks. Although compression of a single eddy 
leads to the first order Fermi acceleration, the first or- 
der effect cancels out after the particle encounters many 
compression and expansions and only the second order 
effect leads to increase of momentum through diffusion 
process (see Bykov & Toptygin 1982). 

In this section we have implicitly assumed that parti- 
cle diffusion time is similar to or larger than the eddy 
turnover time on the outer scale. In the fast diffusion 
limit the our earlier estimates of the turbulent acceler- 
ation are valid. Therefore the the momentum diffusion 
coefficient is of the order ~ p 2 v\/D\\, if the particle does 
not return back to the high density compression. If it re- 
turns before the compression disappears this estimated 
should be multiplied by N > 1. 

6. ACCELERATION BY ELECTRIC FIELD ARISING FROM 
SLOW MODES IN INCOMPRESSIBLE LIMIT 

Let us consider the extreme case of high (3 plasma, 
namely, incompressible plasma. In the incompressible 
limit only two modes, Alfven and slow (pseudo- Alfven) 
exist. Slow modes compress magnetic field. 

6.1. Betatron Acceleration Process 

In incompressible conducting fluid a slow wave causes 
an oscillation occurring in the plane spanned by Bo and 
k vectors, where Bo is the mean magnetic field, or more 
precisely local mean magnetic field, and k is the wave 
vector. Therefore, the electric vector E oc — Vk X Bo is 
perpendicular to both Bo and k. Here, is the ampli- 
tude of the oscillation velocity of the pseudo- Alfven wave. 
The electric field is parallel to the wave front (see Figure 
EJ). Thus, slow modes can create non-zero V x E (see 
Figure H^b)). When a particle gyrates around B , the 
particle feels electric field caused by turbulence. When 
the electric field has non-zero curl, the particle experi- 
ences either acceleration or deceleration. An accelerating 
eddy 7 is shown in Figure Hfc). 

7 An eddy-wave duality in the description of MHD turbulence is 
discussed in CL03. 
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Fig. 1. — Compression rates in strongly magnetized compressible MHD turbulence. {Left) When the sonic Mach number is higher, mass 
concentration is higher. For example, the triangle on the dotted curve, which is for M B ~ 2.2, tells us that about half of the mass occupies 
20% of the volume. The diamond on the dashed curve says that 25% of mass fills 2.7% of the volume. (Right) The compression factor 
may be defined by volume filling factor divided by the mass fraction for triangles and diamonds in the left panel. We plot these values as 
function of the sonic Mach number. Note that the solid line (MHD) and the dashed line (hydro) have similar slope. Although the solid line 
has a shallower slope, this figure roughly suggest that compression rate in isothermal MHD medium shows a dependence on Mach number 
similar to isothermal hydrodynamic one. 
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Fig. 2. — Slow modes and non-zero V X E. (a) Direction of E. Electric vector is perpendicular to both Bo and k, which means that, 
when viewed from the top, electric field is parallel to the wave front, (b) Nonzero V X E by slow modes. The dashed and the dotted lines 
denote wave fronts of wave 1 and 2, respectively. Arrows are the directions of electric field. Two slow waves can create non-zero V X E. 
(c) Acceleration by an elongated slow mode eddy. The particle experiences constant acceleration in the eddy. The integral /(V X E) ■ ds 
should not be zero in the eddy. 



Consider a particle with gyro radius r g moving through 
turbulent eddies of (perpendicular) scale I > r g . Then 
the average electric field the particle feels during one gyro 
rotation is 

E = — !— I E • dl = — !— / (V X E) • ds 

2-Kr g J 2-KTg J 

= t^— [ (B V • v, - B • Vvj) • ds, (44) 
2nr g c J 

= -— !— /(B • Vvj) • ds, (45) 

2-KTgC J 

where dl is the infinitesimal line clement vector, ds the 
infinitesimal area element vector, E = — v x Bo/c, and c 
the speed of light. The area element vector ds is taken to 
be parallel to B . Eq. I|44|) is valid for both incompress- 
ible and compressible fluids. However, eq. I|45|l is valid 
only for incompressible fluids because we used V • v = 



when we derived it. The integral is zero for Alfven modes 
because v; • ds = 0. For slow modes, the integral is ap- 
proximately 

J (B • Vv,) • ds « Bo^H {TTr 2 g )/l\\ , (46) 

where ujii is the velocity component of pseudo- Alfven 
modes parallel to Bo- Thus, we have 

E^-B v l \\r g /2l\\c. (47) 

The amplitude of E is the largest 8 at I ~ r g . Indeed, 
when I > r g , since ujy ~ vi oc I 1 / 3 and Zy cx l 2 ^ 3 (GS95), 

E oc l^ 1 / 3 and E increases with the decrease of /. On 
the contrary, when I < r g , the particle traverses many 

_ 8 Note, that for particle acceleration not only the amplitude of 
E but also its correlation time matters. The latter is larger for 
larger eddies. 
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uncorrelated eddies during one gyro orbit. The value of 
(V x E); in an eddy of size I is of order (V x E); ~ 
Bo • Vv;/c ~ BqVi/Iuc oc I' 1 / 3 . This value times As ~ 
ttI 2 is the total contribution to E by the eddy (c.f. eq. 
H46|) ~). There are ~ {i^/l 2 ) 1 / 2 uncorrelated eddies inside 
the particle's gyro orbit. Therefore the area integration 
(cf. eq. H45[) ~) of these small scale fluctuations over a circle 
of radius r g > I yields E oc (r^/Z 2 ) 1 / 2 ^ 2 )^ 173 
As the result, the E increases with I, 



r g l 2 ' z . 



6.2. Structure of electric field of slow modes 

In this section we use a data cube obtained from a 
direct incompressible MHD turbulence simulation with 
256 3 grid points (Cho & Vishniac 2000) to study the 
structure of E. The data cube clearly shows that MHD 
fluctuations do produce non-zero V x E. 

In Figure Ola), we plot contours of equal second or- 
der structure function for electric field in a local frame, 
which is aligned with the local mean magnetic held Bx: 
SF 2 (r h r ± ) =< |E(x + r) - E(x)| 2 > avg . over x , where 
r = ri|fi| + r±r± and f || and fj_ are unit vectors parallel 
and perpendicular to the local mean field Bl, respec- 
tively (see Cho & Vishniac (2000) for a detailed discus- 
sion of the local frame). The contour plot clearly shows 
the existence of scale-dependent anisotropy: smaller ed- 
dies are more elongated. By analyzing the relation be- 
tween the semi-major axis (~ hi ~ 1/&||) and the semi- 
minor axis (~ l± ~ l/fcj_) of the contours, we can obtain 
the relation between fc|| and k±. Here k± and fc|| are 
wave numbers perpendicular and parallel to the mean 
magnetic field, respectively. The result in Figure HHb) is 

2/3 

consistent with the GS95 type anisotropy: fcii oc . In 
Figure we plot § r E • d\ for r g =8 in grid units. The 
particle orbit lies within the inertial range of turbulence, 
which spans from ^100 grid units to a few grid units 
in the simulation data. We can see that structures are 
elongated along the mean field direction. 

6.3. Adiabatic gain/loss in parallel directions 
The rate of momentum change by the electric field is 



dp± 
dt 



qE 



qBpvgrg p±cvi\\r g p±vg 



2/n 



(48) 



where we used qBo = p±c/r g , where p± is the particle 
momentum perpendicular to magnetic field. Since the 
integrand in eq. is B Vi_ • v;j_, where '_L' denotes 
the component perpendicular to Bo, we can rewrite eq. 
(gSJ) as follows: 

dpx/dt = -piV± ■ v_l/2. (49) 

This means that the acceleration by electric field is 
achieved when large scale turbulence compresses mag- 
netic field in perpendicular directions. The change of 
perpendicular component of energy is 



El 

Ef 



9 P± _ 

_ 7 ' 



(50) 



where E is the kinetic energy and we consider non- 
relativistic case for simplicity. 

When we consider incompressible fluids, compression 
in perpendicular directions causes expansion in parallel 



directions, which in turn results in adiabatic energy loss. 
The rate of momentum change by the parallel motions 
alone is 

p\\ ~ h ' 



(51) 



Therefore, the change of parallel component of the energy 
is 

E« 



E'< 



'II p W 

From equations i|50|) and 152fl . we get 



2^ = -2 V f. (52) 



d{E£+E() 
dt 



K, V l\\ 



{Ef-2Ef)^ = Q 



(53) 



when particle distribution is isotropic (i.e. E^_ = 2Efr). 

Therefore, if the particle distribution is isotropic, the 
second term in eq. (|44fl . or the term in eq. I|45|). is 
canceled by adiabatic gain/loss in parallel directions and 
only the V • v ; term in eq. survives, the result of 

which is equivalent to the hydrodynamic approach of P88 
that we discuss in the next section, i.e. 



dp/dt ~ — pV • v. 



(54) 



However, it is clear from eq. (|53|l that an acceleration 
instability is present when betatron acceleration acts. In- 
deed, as electric field increases E±, the gain of the en- 
ergy gets positive. As processes of scattering tend to 
isotropize particle velocities, we may assume that this 
instability acts on the time scale of randomization of a 
particle. This time can be of the order of mean free path 
divided by the particle velocity and multiplied by a fac- 
tor larger than unity (see discussion in Jokipii 1968). For 
our simplified treatment we conservatively assume that 
this factor is order of unity. In other words we claim 
that, on the mean free path scale of the particles, non- 
resonant effects are important for particle motions and, 
therefore, the hydrodynamic approach of P88 fails. 

Consider an eddy a few times smaller than the mean 
free path. Let this scale be V . Suppose that the eddy 
is compressed in the perpendicular directions when a 
particle is back-scattered in the eddy. During the back- 
scattering, the particle's pitch angle is obviously large: 
around ~ 90°. Suppose that the particle stays in this 
large pitch angle state for At. Later, the particle will 
have smaller pitch angles. On average, the particle will 
be in the small pitch angle state for At. While the 
particle is in the small pitch angle state, the particle 
will travel more distance because the parallel velocity is 
higher. Therefore, the particle will pass through several 
eddies of the scale I', the acceleration/deceleration by 
such eddies will roughly cancel out. As a result, the mo- 
tions of eddies that contain the turning points are more 
important. Since such motions are mutually indepen- 
dent, we can have diffusion of momentum. 

Since large scale eddies do not provide net acceleration 
when particle distribution is isotropic, we do not con- 
sider large scale eddies in the incompressible limit. In 
the next subsection, we consider acceleration by eddies 
on the mean free path scale. 

6.4. Estimates of D p 
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When particle distribution is isotropic, we do not have l|53(l is no longer valid on the mean free path scale, 
net acceleration (see eq. IpBl) by large scale eddies. How- Here, we use properties of strong MHD turbulence, the 

ever, as we discussed in the previous subsection, equation simplest case of which is that turbulence is isotropic on 
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the energy injection scale L and that Bq ~ 6^, where bi, 
is the rms fluctuating field. For eddies whose parallel size 
is comparable to the mean free path by pitch-angle scat- 
tering, l m f p , particles are accelerated for At ~ l m fp/ v pth 
where v p a is particle velocity. Note that in the absence of 
fast modes the mean free path may be very large as scat- 
tering by Alfven and slow modes is marginal (Chandran 
2000; YL02). The momentum change over the particle's 
eddy crossing time, At ~ l^/v pt i, is 



Ap < 



(55) 



(dp±/dt)At ~ p±v l \\/2v p ti, 

where we used the acceleration that arises from magnetic 
field compression (see eq. (@H10 : dp/dt ~ pvi^iow/U (see 
Appendix). Therefore, 



jyincoinp 
V 



(Apf/At ~ ( P 2 vl slow /6v ptl L), 

^ p 2 V A ^ VL,slo 



L 



VA 



(56) 



where we used vq ~ vl,sIowQ I 'L) 1 /^ and l\\ ~ L 1 / 3 / 2 / 3 . 
For vl.si ow ~ u^, this is ~ va/v p u times smaller than 
eq. (30) in C03 (see Appendix A). The ratio va/v p u 
is small for the interstellar medium, since va is a few 
km/sec and v pt i can be as large as the speed of light. 
However, for a typical solar coronal loop, this ratio is ~ 
1/30 since va ~ 10 4 km/sec (see C03). This mechanism 
is more efficient than particle acceleration by large scale 
compressive motions by slow modes (§2J) in the high-/? 
limit, where the efficiency of the latter goes to zero. 

As the interaction happens over the time of ballistic 
passage of a particle through an eddy it is fast for all the 
physically interesting cases. Random walk in terms of the 
applied field may also be applicable, if turbulence scale 
is less than the mean free path of a particle. This would 
only decrease the efficiency of the mechanism, however. 

7. ACCELERATION BY WEAK TURBULENCE 

Although we believe that for most of astrophysical, e.g. 
interstellar (see Armstrong, Rickett & Spangler 1995), 
conditions the GS95 scaling of slow modes presents the 
best-known fit, weak turbulence (Galtier et al. 2000) 
may also arise in astrophysical situations (Saur et al. 
2002). In weak turbulence, vl <C va and turbulent dif- 
fusion is very slow. Particle acceleration by anisotropic 
weak turbulence is useful, as MHD perturbations at large 
scales may evolve initially along weak turbulent cascade 
before turning into the strong cascade (see a discussion 
in Cho, Lazarian & Vishniac 2003). In this section, we 
consider particle acceleration in weak turbulence. Since 
large-scale Alfvenic motions do not accelerate particles, 
we discuss slow and fast modes. We assume weak turbu- 
lence in compressible medium follows similar scalings as 
in incompressible one, which is a non-trivial conjecture 
for low-/3 medium. 

7.1. Incompressible limit: betatron acceleration 

In incompressible limit the slow modes become pseudo- 
Alfvenic perturbations that move with the Alfven ve- 
locity. We should keep in mind that for weak tur- 
bulence slow mode perturbations have velocities sub- 
stantially less than the Alfven velocity. To calculate 
D p in weak turbulence, we can still use the relation 
Ap ~ (dp/dt) At — (pvi,siow /l\\)(l\\/ Vpti) and At — l\\/v pt i 



(see eq. Q55[l). However, anisotropy and scaling of vi\\ are 
different. Weak turbulence has an extreme anisotropy, 
Z|| = constant, and v t ~ I 1 / 2 (Galtier et al. 2000). We 
have 



j^incom,weak ^ (Ap) 2 /At 



P v L,slow{ l mfp/ X)/(6v pt ;L||) 
%fpVA 



p 2 V A ( V L ,slow \ 2 l r 



VA 



(57) 



where we used Zii 



L 



L\\Vpti 
This is smaller than 



D incomp,strong m eq Q by a factor f _ l mfp /L if 

L||~L. 

7.2. Slow modes: fast diffusion limit in compressible 

fluid 

The acceleration processes by weak and strong turbu- 
lence in the fast diffusion limit are very similar. As in 
§4.1 consider Dm Lc s , where c s is the propagation 
speed of the slow wave and depends on the plasma (3: 
c s ~ va for high-/? and c s ~ a (=sound speed) for low- 
/?, where we ignored the angle dependence for simplicity. 
Inserting 

Ap ~ (pv liS i ow /l\\)(l 2 /D\\) ~ pvi fS i ow l\\/D\\ (58) 



and 



At ~ q/Du 



into D p ~ (Ap) 2 / 'At, we get 



d p,i ~P v Uslow /D h 



(59) 



(60) 



which is largest on the outer scale. 

When we take into account multiple re-entry of the par- 
ticle to the same eddy, we need to multiply D p above by 
N w ~ ram{^ D\\Ti,rand/li Minima), where n,rand is the 
eddy randomization time of the slow modes on scale I. In 

weak turbulence, n,rand ~ X^hwave ~ ^ed&ylKviave. oc 

I, where x = bil\\/B a l ~ ti, wave /ti, e ddy (GS95) and wc 
used U^ddy ~ l/vi ~ l 1/2 and t^ wave = L/c s = con- 
stant. Here bi is strength of the fluctuating field on scale 
Therefore the first term in the parenthesis is propor- 
tional to ~ Z -1 / 2 . On the other hand, the second term, 
witch describes the characteristic scale of field line di- 
vergence, should decrease when I decreases. The reason 
is as follows: the separation of field lines will show a 
random- walk behavior on scales larger than L. For a sin- 
gle eddy on the outer scale, separation of two adjacent 
field lines is ~ Lub/Bg, which will be the average path 
of the random walk. To reach the distance of L, we need 
~ (LB /L\\b) 2 ~ Xl 2 steps, where \h is X evaluated on 
the outer scale. Therefore, the second term in the paren- 
thesis, Mi^ s i ow , evaluated on the outer scale is ~ Xj} ■ 
On the other hand, weak turbulence will show transition 
to strong turbulence deep down in the inertial range. We 
know that the factor Mi >s i ow is of order unity when tur- 
bulence is strong. To summarize, the factor Mi tS i ow is 
very large on the outer scale and order unity deep down 
in the inertial range. The most logical conclusion from 
this is that the factor Mj jS ; om decreases as I decreases. 
As the result, the product of the momentum diffusion 
coefficient in eq. ||dT7|) and the factor N w decreases as I 
decreases. 
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Therefore, the efficiency of particle acceleration is 
largest on the outer scale. On the outer scale, the factor 
N is 



D \\X L 2t L,wave/L,X L 2 } 
~ minjXi 1 \Jt L ,wave/tL,diff^Xl 2 } (61) 

Note that \l <C 1 in weak turbulence. The first term 



is smaller when tL,diff > X 2 *£,' 
coefficient is 



The net diffusion 



D p ~P 2 {v 2 LMow /D\{)N w , 



(62) 



where the factor N w ^> 1. 

All in all, in the fast diffusion limit the acceleration 
efficiency is increased by the factor of JV„ which is larger 
than the factor N for the strong turbulence. However, 
while for strong turbulence the injection velocity may 
be of the order of Alfven one, for weak turbulence the 
injection velocity should be smaller than va- 

7.3. Slow modes: slow diffusion limit in compressible 

fluid 



In this subsection, we consider Z)|| <C Lc s 



As we dis- 
cussed in earlier sections, as long as Du l\\c s , random 
walk-like behavior is suppressed. The random walk-like 
behavior is fully recovered only for the eddies that satisfy 
Z?|| ~ ii|C s . For slow and fast modes in strong turbulence, 
it is possible to find such eddies. However, since Z|| is con- 
stant in weak turbulence, it is not possible to find such 



a scale that satisfy D 



|C a . 



Therefore, acceleration 



of particles by large-scale motions is inefficient in weak 
turbulence when diffusion is slow. 

As we mentioned earlier, down the inertial range weak 
turbulence undergoes a transition to the strong one as k± 
increases. The acceleration by such turbulent motions is 
given by eq. Q25J1 f° r low-/3 plasma and eq. (|28Jl for 
high-/3 plasma where the velocity of slow modes at the 
transition to the regime of strong turbulence should be 
used. 

7.4. Fast modes: slow and fast diffusion limits 

Acceleration that is induced by fast modes in weak 
turbulence is analogous to the acceleration by fast modes 
in strong turbulence. Therefore the momentum diffusion 
coefficient in fast diffusion limit is given again by eq. 
I|31|) . but instead of multiplying it by the factor N, as in 
the case of strong turbulence we should multiply it by a 
larger factor N w to account for a higher probability of a 
particle to return back to the compressive eddy. For slow 
diffusion limit, the contribution by large scales eddies is 
given by eq. 1)1 0|l and therefore sub-dominant. This is 
due to the fact that for weak turbulence the turbulent 
diffusion is slow. The actual contribution to acceleration 
in the slow diffusion limit is made by eddies with size 
I ~ D\\/cf. The amplitude of the fast modes depends on 
the properties of driving. 

8. ACCELERATION BY PITCH-ANGLE SCATTERING 

Here we consider acceleration by pitch-angle scatter- 
ing. We do not need to know detailed physics of pitch- 
angle scattering. However, we assume that characteristic 
propagation speed of the scattering agent is ~ va ■ When 
the other velocities are involved, appropriate correction 



is needed. Pitch-angle scattering is possible by both com- 
pressible and incompressible motions. In most cases, this 
is the scattering that determines the particle mean free 
path and the Du . 

When a particle interacts resonantly with small-scale 
waves, both the pitch-angle, the angle between momen- 
tum p and Bo, and the particle momentum change 
(Jokipii 1966; Skilling 1975; Schlickeiser & Achatz, 1993; 
see also Shalchi et al. 2004). This pitch-angle scattering 
results in a diffusion process in momentum space and the 
particle distribution function / follows 



df d_f D 0/ 



dt 



JL( D 21 

d/z \ w dp 



p 



l__d 

2 dp 



P 



df 



dl 

pp d P 



(63) 



(Schlickeiser & Achatz 1993), where fj, — cos 9 and 9 is 
the pitch-angle. 

The standard assumption for the studies of the accel- 
eration by large scale compressions (see P88, C03) is that 
the diffusivity arises from the the pitch-angle scattering, 
while the D pp arises from the large-scale compressions. 
However, it is easy to see that this is not true. Indeed, 
the pitch-angle scattering frequency is 



v = 2£> W1 /(1 - M 2 )- 



(64) 



The /^-averaged momentum scattering coefficient D p is 



D p 

p 



itch 



d[iD p 



(65) 



for isotropic particle distribution. In general, the mo- 
mentum diffusion coefficient D pp is ~ p 2 {vA/vpti) 2 
times the pitch-angle diffusion coefficient D^: D pp ~ 
p 2 (v A /v ptl ) 2 D^ (Skilling 1975; Chandran 2000; YL02). 
In fact, there is a simple interpretation for this: during 
each back-scattering, the fractional momentum change 
Ap/p is ~ va/v p u. Therefore, roughly speaking, 



Df ch ~ p 2 (v A /v p uY 



dfiD^ ~ p 2 {v A /v pt i) 2 v, 



(66) 

where we assumed J dfiD^^ ~ J dfiD^^/(l — /i 2 ). More 
precisely speaking, 



j-^pitch 
P 



Ap 2 



VA 
Vptl 



(67) 



where the angled bracket (...) indicates an average over 
the pitch-angle cosine /i and v + and ^_ are the pitch- 
angle scattering rates by waves moving parallel and anti- 
parallel to Bo, respectively (Skilling 1975). This equa- 
tion reduces to eq. I|66l) when u + = which means 
that parallel and anti-parallel waves have equal power. 
In other words, eq. (|66|l is valid when turbulence is bal- 
anced. Note that D p is smaller than that in eq. I|66|l 



when turbulence is imbalanced. Noting v 



Vpti/lmfp, we can rewrite eq. (16611 as 

Df ch ~p 2 {v A /v ptl ?v 

~ P 2 {VA/Vptl) 2 v p tl/I m fp ~ P 2 V 2 A / (Vptllmfp) 

^ P 2 V A Lv A 

L ImfpVptl 



4ti/ D \\ 



(68) 
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Alfven, slow and fast modes contribute to pitch angle 
scattering. However, due to inefficiency of scattering by 
Alfven and slow modes the corresponding l m f p is very 
large. Therefore the fast modes dominate pitch-angle 
scattering. The pitch-angle scattering is comparable to 
slow or fast modes acceleration by large scale compres- 
sions in fast diffusion limit and dominates them in slow 
diffusion limit in the absence of shock compression. How- 
ever, in case particles are scattered by a static magnetic 
field, e.g. magnetic mirrors created by molecular clouds 
(see Chandran 2001), the pitch-angle scattering will not 
accelerate particles efficiently. 

Pitch-angle scattering also occurs in weak turbulence. 
To calculate D p , we can use eq. However, the 

pitch- angle scattering efficiency, hence Imfpi is l ess than 
that in strong turbulence because anisotropy is more pro- 
nounced and fluctuation of magnetic field is smaller in 
weak turbulence. Weak turbulence in fast diffusion limit 
is dominated by large-scale non-resonance acceleration 
given by eq. (ffi2"|) . 

Transit Time Damping (TTD) and gyroresonance both 
produce pitch angle scattering 9 . 

Gyroresonance is based on the resonance between the 
frequency of gyrorotations and the frequency at which 
the particle meets bumps in the magnetic field. There- 
fore sufficiently small scales matter. The TTD is essen- 
tially Cherenkov-type interaction, which allows interac- 
tions of particles with large-scale turbulence. As the re- 
sult, TTD is important for the situations when turbu- 
lence is damped at small scales. Note, that TTD, unlike 
gyroresonance, requires compressions and therefore is in- 
efficient for Alfvenic perturbations. 

At any rate, we expect that the effect of TTD is very 
similar to that of gyroresonance scattering 10 . For in- 
stance, D pp is ~ p 2 (v A /v p ti) 2 times the pitch-angle dif- 
fusion coefficient D^: £>J p TD ~ p 2 (v A /v pt i) 2 D^ D (Yan 
& Lazarian 2004). 

A frequent argument against the TTD is that such 
mechanism accelerates the particles along the magnetic 
field and therefore gets quickly saturated as the fraction 
of modes available for acceleration decreases. Two effects 
act to counteract this, however. First of all, magnetic 
field wandering (see Lazarian & Vishniac 1999) changes 
the direction of magnetic field in respect of wavevectors of 
fast modes. Moreover, the anisotropic distribution with 
particles moving mostly along magnetic field lines is sub- 
jected to instabilities that will randomize particles. We 

9 Formally speaking both TTD and gyroresonance are given by 
the same condition (see Melrose 1968): u> — kuVn — nQ/7, where 
frequency w and a wavevector component parallel to magnetic field 
fen are the characteristics of the magnetic turbulence, while and 
Q are, respectively, the component of particle velocity parallel to 
magnetic field and gyroresonance frequency of the particle. The 
peculiar feature of the TTD is that it corresponds to n = and 
therefore does not depend on particle frequency. This entails sub- 
stantial differences. In terms of calculations, for instance, while 
the contributions from large-scale MHD motions are unphysical for 
gyroresonance, as those would violate the adiabatic invariant con- 
servation (Chandran 2000; Yan & Lazarian 2003; Shalchi, Yan & 
Lazarian 2005), they should be accounted for in TTD calculations 
(YL02). 

10 TTD, unlike gyroresonance, does not randomize the distribu- 
tion of particles. Instead, it increases the parallel particle velocity 
and does not change the perpendicular velocity. However, this 
eventually should result in instabilities that randomize particle ve- 
locities. 



note that, while TTD may be also driven by large-scale 
compressions, it is different from the large-scale stochas- 
tic mechanisms we deal for the most part of the paper. 
It does require a rather coherent interaction with the 
Fourier component of the random field. A particle surfs 
over the wave rather than bounces back and forth within 
a compressed region. 

9. COMPARISON OF ACCELERATION MECHANISMS 

Above we discussed different mechanisms of particle 
acceleration that make use of turbulence. These mech- 
anisms can be roughly divided into two groups: pitch- 
angle, which include gyroresonance as well as TTD and 
non-resonant mechanisms that rely on large-scale com- 
pressions. Dealing with the second group of mechanisms 
we obtained the acceleration rates for compressible (fast 
and slow) modes in both weak and strong turbulence, 
considered limits of slow and fast particle diffusion. In 
addition, we showed that, the compressions of magnetic 
field in incompressible fluid can still provide acceleration 
of particles over their mean free path. 

We summarize the results for strong MHD turbulence 
in Tabic 1. When v A ~ v L ^ ow ~ v L j ast , we have 

D vncom < pslou, < pfast < ppitch _ (gg) 

Relative magnitude of D^ tch and rj^ u P ersomc varies. In 
slow diffusion limit (i.e. t^^iff > th.eddy), depending on 
the value of Qtd, D p dow can be larger than D^ ast or 
D pitch_ However, since the factor Qtd is no larger than 
a logarithmic factor, we do not consider the possibility 
seriously. Noting 

(_^_) „ (tj^ll) , (70) 

\l-mfpVptl J \tL,eddy J 

we find that the momentum diffusion coefficient is no less 
than 

D p ~ p 2 v A /L (71) 
for slow diffusion limit. In this limit, jyi£ tch is larger than 
D siow Qr D f p ^. On the other hand, in fast diffusion limit, 
D s p low , Df p as \ and DP ltch are all similar and they can be 
less than p 2 v A /L. 

In slow diffusion limit, when we ignore the Qtd fac- 
tor, slow modes are less efficient than fast modes. This 
is because slow modes have longer wave periods (see Ap- 
pendix B). As we mentioned earlier, D p ~ p 2 (Ap) 2 / At. 
Let us consider j3 ~ 1 case. For slow modes Ap ~ 
p(V ■ vt)At ~ p(vi/l\\)(l\\/v A ) ~ pvi,siow/v A . For fast 
modes Ap ~ p(V ■ vj)Af ~ p(vi/l)(l/v A ) ~ pvij ast /v A , 
where we used c/ ~ v A . Therefore, there is no big differ- 
ence in (Ap) s i ow and (Ap)f ast . However, they have dif- 
ferent At. For slow modes, we have At ~ l\\/v A (GS95). 
For fast modes, we have At ~ l/v A because fast modes 
are wave-like. Since lu 3> I on small scales, fast modes 
have much shorter At on small scales. Therefore, fast 
modes are more efficient. 

Since pitch-angle scattering is comparable to or more 
efficient than slow or fast modes, we estimate the value of 
jjpitch f Qr var i ous astrophysical objects. We first rewrite 
the momentum diffusion coefficient by pitch-angle scat- 
tering as follows: 

D p it cH „ P^A f tLMt \ (?2) 
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where tL, W ave ~ L/va- From the relation D p ~ 
(Ap) 2 /At <~ p 2 /r acc , we can define the acceleration time 
as follows: 

Tacc = "p- (73) 

which is the timescale for doubling the momentum. The 
results are given in Table 2. 

10. DISCUSSION 

Our calculations in the paper above reveal that in fast 
diffusion limit, namely, when a particle diffuses away 
from a compression on the time scale less than the rele- 
vant crossing time, weak and slow turbulence accelerates 
particles in very similar way. The most notable difference 
between the two types of turbulence is that the particle is 
more likely to return back to the compression created by 
weak turbulence. At the same time, in the slow diffusion 
limit the differences in acceleration by weak and strong 
turbulence are substantial. They stem from both the dif- 
ference in turbulent scaling and from the differences in 
turbulence diffusivity. 

Our study allows to answer a number of important 
questions, e.g. we can answer the question we posed at 
the beginning of our paper. For instance, similar to the 
case of scattering discussed in YL02, fast modes should 
dominate particle acceleration (cf. C03). We find that 
in many cases small scale turbulent motions (comparable 
with the particle gyroradius) that scatter particles also 
accelerate the particles more efficiently than the turbu- 
lent motions at large scales. TTD which is also efficient 
for acceleration of particles but utilizes both small scale 
and large scale compressions does not explicitly rely on 
particle randomization through scattering. Unlike the 
TTD, the stochastic mechanisms that deal with large 
scale motions do rely on scattering. Therefore, in many 
instances it is enough to limit the studies of MHD tur- 
bulent acceleration to gyroresonance and TTD (see Pet- 
rosian, Yan & Lazarian 2005). The particular cases when 
alternative mechanisms of stochastic acceleration are im- 
portant will be discussed elsewhere. 

Similarly to the study in Yan & Lazarian (2004) we as- 
sumed above that the turbulent energy is being injected 
at the large scales and it cascades to small scales. One of 
us (AL) studies currently a possibility of an additional 
energy injection that happens at small scales and the 
effects of this on cosmic ray. For instance, streaming in- 
stability (see Cesarsky 1980) may be a source of injection 
of energy in the form of waves (see YL02, Farmer & Gol- 
dreich 2004, Yan & Lazarian 2004). In this case, Alfvenic 
motions would be more efficient in scattering compared 
to the predictions in YL02. Moreover, the gyroresonance 
acceleration by such motions will be also efficient. 

For the scattering calculations in Yan & Lazarian 
(2004) the plasma damping effects were systematically 
taken into account and proven to be very important. In 
our present paper the damping of turbulence is only im- 
plicitly present through our defining of the scales where 
turbulence is present and through the mean free path of 
the particle l m fp- For many of the discussed mechanisms 
the damping is unimportant as they are based on the in- 
teractions with large-scale fluctuations. However, for the 
gyroresonance the damping is important and the results 
from Yan & Lazarian (2004) can be applied directly. 



The results that we obtained in the paper are order 
of magnitude only. Although calculations for gyrores- 
onance and TTD acceleration using quasi-linear theory 
are available (see Yan & Lazarian 2004) we did not use 
more precise formulae. The merit of our work is that 
it provides a rough comparison between different mech- 
anisms of stochastic acceleration, allows us to define the 
dominant processes of stochastic acceleration, and deter- 
mines the necessary conditions for the different mecha- 
nisms to operate. We believe our rough treatment may 
be justified and considered adequate for many purposes 
in view of various uncertainties of the model of turbu- 
lence. For instance, an uncertainty arises from determin- 
ing the amplitudes of turbulent perturbations. Indeed, 
the acceleration is caused by slow and fast modes, while 
usually only total velocity dispersion is available through 
observations (see Lazarian & Pogosyan 2000, 2004 and 
references therein). Yet extra uncertainties stem from 
the uncertainty of the the injection scale and the scale of 
transition from weak to strong turbulence, if this transi- 
tion takes place. 

We hope that some of the uncertainties can be clari- 
fied. For instance, potentially compressible and incom- 
pressible motions can be separated by statistical analysis 
of spectral line variations (Lazarian & Esquivel 2003, Es- 
quivel & Lazarian 2005). Observational studies of mag- 
netic field intensities and further progress in theory of 
MHD turbulence should make the transition from weak 
to strong turbulence more certain. Therefore a more pre- 
cise treatment of the acceleration processes would be wel- 
come in future. 

In the paper we discussed the acceleration by Alfven, 
slow and fast MHD modes. The possibility of considering 
these MHD modes separately from the rest of the MHD 
cascade stems from the scaling formulae in CL03 that 
shows that the energy exchange between different MHD 
modes drops along the cascade. 

Our results show that some issues in the earlier studies 
of turbulent acceleration require revisions. For instance, 
the non-resonant acceleration in the slow diffusion limit 
is in general more efficient than in P88. In addition, 
fast modes are usually more important for acceleration 
than slow modes (cf. C03). Moreover, it is not good 
to disregard resonant acceleration while considering the 
non-resonant one. The scattering that is required by the 
non-resonant acceleration provides the acceleration rates 
that are usually comparable or larger than those arising 
from the non-resonant processes. 

We believe that the importance of the stochastic ac- 
celeration is underestimated in the literature. It looks 
that this sort of acceleration is important at least for 
Solar flares, clusters of galaxies and gamma-ray bursts 
(see Petrosian 2001; Lazarian et al. 2002; Brunetti et 
al. 2004; Petrosian, Yan & Lazarian 2005; Brunetti & 
Blasi 2005). In addition, stochastic acceleration can also 
be important for acceleration of the particles within the 
turbulent reconnection regions (see Lazarian 2005). 

We also mention that the particle acceleration mech- 
anisms discussed above are applicable not only to ener- 
getic particles, but also to charged dust grains. Resonant 
scattering and TTD were applied to charged dust grains 
in Yan & Lazarian (2003). Because of the relatively small 
velocities of dust grains, the factor {va/v p u) 2 is usually 
larger than 1. Therefore the gyroresonance and TTD 
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acceleration are more efficient than scattering and ran- 
domization of grains, and the resonant processes domi- 
nate the MHD acceleration of dust grains. In many cases, 
due to inefficient magnetic scattering, the mean free path 
of the grain is determined by gaseous or plasma damp- 
ing. Such a motion is mostly ballistic with the adiabatic 
invariant of the grain conserved. As a result, the large 
scale compressions change grain velocities by a factor of 
order unity. 

All the results above assume turbulence where random 
motions are homogeneously distributed in space. Turbu- 
lence intermittency should affect the distribution func- 
tion of the accelerated particles. This issue, however, is 
beyond the scope of the present paper. 

11. SUMMARY 

We have studied particle acceleration by weak and 
strong MHD turbulence. We considered both case when 
the particle diffusion is slow and fast compared to the 
life-time of a compressible perturbation. 

1 . Efficiency of particle acceleration depends upon tur- 
bulence being strong or weak at the injection scale 

2. Pitch-angle scattering in most cases dominates cos- 
mic ray acceleration when particle diffusion is slow. 



3. Pitch-angle scattering, slow modes, and fast modes 
give similar acceleration efficiency when particle diffusion 
is fast. 

4. Turbulent diffusion arising from hydromagnetic 
motions allows more efficient non-resonant acceleration 
when particle diffusion is slow. 

5. For strong turbulence, fast modes dominate the 
non-resonant acceleration. 

6. High compressions arising in super-sonic MHD tur- 
bulence increase efficiency of non-resonant particles ac- 
celeration. 

7. Incompressible MHD turbulence can provide non- 
resonant acceleration of particles over mean free paths. 
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APPENDIX 
CALCULATION OF D P 

Dp from anisotropic slow modes 

Chandran & Maron (2004b) considered particle acceleration by strong turbulence. They assumed that the eddy 
turnover time t e ddy ~ L/vl) is less than the particle diffusion time tdiff ~ L 2 /D, i.e. that the particles can diffuse 
out of an eddy before it is randomized. The corresponding At = tdiff and V • V/, ~ Vl/L and obtained (see i j2.1(l 

D p ~ p 2 (v 2 L /L 2 )(L 2 /D) ~ p 2 v 2 L /D ~ p 2 v 2 A /D, (Al) 

where vl ~ va is assumed. They observed that, when particles diffuse along magnetic field lines, they can reenter the 
same outer scale eddy 

N ~ min{yjD\\T rmd /L,M} (A2) 

times 11 before the outer scale is randomized. Here T ran d is the eddy randomization time on the outer scale of turbulence. 
The factor M is defined as z s /L, where z s is the distance along magnetic field lines over which two adjacent magnetic 
field lines get separated by the distance L. The estimates for this factor range from 1 (Narayan & Medvedev 2001) to 
5 or 7 (Chandran & Maron 2004b). Therefore the value of N is no more than a few in strong turbulence. We observe 
that intrinsic diffusion of particles due to scattering should decrease M, while if turbulence is weak, M should increase. 

On the other hand, when tdiff > t e ddy, particles are confined within the eddies until the eddies are randomized. C03 
used V • vl ~ vl/ L and At <~ t e ddy (~ L/vl) for outer scale eddies and obtained 

D P ~ p 2 (v 2 L /L 2 )(L/v L ) ~ p 2 v L /L ~ p 2 v A /L (A3) 

(see also H2.1fl . C03 claimed that the small scale structure of slow mode turbulence is also important when d m i n VA *C 
D\\ <C Lva, where va is the Alfven speed and d m i n is either the parallel size of eddies at the dissipation scale or 
the mean free path of the cosmic rays, whichever is larger. Particle diffusion time is smaller than the eddy turnover 
time at the outer scale L and larger than that at the scale d m i n . For eddies whose parallel size is larger than 
D\\/va, particles are confined within the eddies until the eddies are randomized. For such eddies, we note that 

V • v/. s i oto ~ (vi lS iow/l)(i ■ £s) ~ vi t aiow/l\\ (see Appendix Bl for details) and At ~ 1\\/va (GS95). Therefore (see eq. 

&), 

D™ 3 ~ (Ap) 2 /At ~ ( P 2 v 2 slow /vAh) ~ (p 2 < s Io W /vaL), (A4) 

where we used vn\ ~ Vl, s Iow(1/ 'L) 1 ' 3 and hi ~ L 1 / 3 / 2 / 3 . Each eddy whose parallel size is between Du/va and L makes 
an equal contribution to D p . Since there are roughly J dl/l ~ In (Lva/Du) eddies in the parallel size range, the sum 

11 The first term on the right is due to diffusion process. For a diffusion process, we can use the relation, d? ~ DyAt, where d is the 
average distance the particle travels through diffusion. During one eddy turnover time of the outer scale, L/vl, the test particle can travel 

d ~ W DuL/vl, . The number we are looking for is d/L ~ • / D^/Lvl ~ J DuL/vl/L. 
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TABLE 1 

Summary of D p 



Dp Estimates Properties of fluid Eq. 

incompressible fluid 1561 

slow modes, low-/? a ' , tL ; diff ^ ^L,eddy 

slow modes, low-/3, t L>di ff > t Lyeddy J2£j 
slow modes, high-/?, t Lidiff < t L ^ ddy 
slow modes, high-/?, t L}diff > t L>eddy OH) 

fast modes, all-/3, t Ltdi ff < t LiWave £HJ 

fast modes, low-/3, th,diff > tL,wave 

fast modes, high-/3, t L}diff > t L>n]ave 

pitch-angle scattering 1681 

super-sonic fluids d ) 



a 'The plasma j3 is defined by the ratio of the gas pressure and the magnetic pressure: (3 = Pg/Pmag (= o? ' /v 2 A 
for isothermal gas), where a is the sound speed and va the Alfven speed. The value of rrt is related to the 
scaling of fast mode: v lJast ~ v L j ast (l/L) m . 

b 'Eddy turnover time at the outer scale L is t^ ^dy ~ ^l v A\ diffusion time at the scale is tj^ di jj ~ L 2 /D|| ~ 
L 2 /ImfpVptl ; wave period at the outer scale is ti, i( iiff ~ ^/ c /> where Cf is the propagation speed of fast modes. 
Note that t L:di f f /t Lt< , ddy ~ iii^/Dy ~ Lv A /l m fpV ptl ; t Ltdif ; /t LtWave ~ Lcf/D^ ~ Lcf/l mf pVp U . Here Z m/p 
is the mean free path of particle back-scattering by pitch-angle scattering. 
c ^Qtd ~ 1 when turbulent diffusion is slow and Qtd ~ h^Z/i^/Dy) when it is fast. 

d 'The parameter 7 is related to the compression rate of the fluid: M] is the ratio of the density in a typical 
compressed region and average density, whereM s is the sonic Mach number. We expect 7 < 2 in MHD 
supersonic turbulence, while 7 ~ 2 in hydrodynamic counterpart. 
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TABLE 2 

Estimates of D„ 



Astrophysical fluids 



Estimates of Df tch a > 



Estimates of T aC c b ^ 



Solar flares 
ICM 

Gamma-ray bursts (30 s~ 1 )p 2 ^ 
Galactic halo c ) 



( 2 Op 2 (t5^7^)( 



5000fcm \ / t L,diff 



°' 5 S ( J ) { 



(3 x ID"* y»"V {-^^) (») 3 x 10? yrs (_^_) (ta) 

lOOOfem 



5000fcm 7 V t L,diff , 

t L ,wave 



3xl0°fcm/sec y 

(10-6 yrs-i)p2 (-^^ 



lOOpc ' 



' t L,diff s 
' t L,diff ' 



( L ) 




\ lOOOfem y 


V tL.diff 1 



1 100km/sec\ 


( L \ 




\ V A J 


\ toopc y 


V t L,diSl J 



a) n pitch p 2 v A I t L,diff ) 

P L \tl,wave) 

b) T = 

»« — -Dp 

c ^ For Galactic halo, / m / p ~ 10 20 cm (Yan Sz Lazarian 2004) for particles with energy larger than ~ lGeV. The diffusion time is 
~ L 2 1 (l m fpc) ~ (100pc) 2 /(/ m /pC) ~ 1000 yrs. The wave period is ~ L/va ~ lOOpc/100/cm/sec ~ 10 s yrs. Cosmic rays in Galactic 
halo arc in fast diffusion limit. We have tL fWa p & /tL.diff ~ 10 3 and r acc ~ 10 9 yrs. 



of all contributions from such eddies is 

-.C03 P 2y A f VL.slow \ 2 i n ( LVA^ 

6L \ va ) n V D\\ 

which for t>L,sZotu ~ V A is equal to eq. (30) in C03. Here we ignored constants of order unity. Note that D^' 03 given 



D C03^i^A(^£^) ln(^), (A5) 
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above does not go to zero when D goes to zero. Therefore, this result is in disagreement with Ptuskin's result in eq. 
©■ 

Acceleration by electric field arising from slow modes in incompressible limit 

In we showed that there is an exact cancellation between energy gain/loss by electric field in the perpendicular 
directions and adiabatic gain/loss in the parallel directions when particle distribution is isotropic. Suppose that, 
somehow, particle distribution is anisotropic. In this case, we can use eq. Ij48(l for dp/dt and 1\\/va for Ai. Therefore, 
the resulting D p will be very similar to those in fJ21 which means that particle acceleration will be efficient if particle 
distribution is anisotropic. 

SCALING OF SLOW AND FAST MODES 
Calculation of V • v 

In Fourier space, V • v becomes ik- Vk£, where i 2 = —1, k is the wave-vector, Vk is the Fourier component of velocity, 
and £ is the unit vector along the direction of Vk- Cho & Lazarian (CL03) showed that 

is oc (-1 + a - /D)A;||k|| + (1 + a - VD)k±k±, (Bl) 
if oc (-1 + a + VD)k\\i\\ + (1 + a + </D)k±k±, (B2) 
where a — a 2 jv\ = /?/2 for isothermal gas, D = (1+a) 2 — 4a cos 2 6, cos 9 — feii/fc, and k± is wave number perpendicular 
to the mean field and ku parallel to it. Hatted vectors are unit vectors. The slow basis i s lies between k|| and —9. 

The fast basis if lies betw een k ± an d k. Here overall sign of £ s and if is not important. 
When /? — > 0, equations (|B2(I and (|B1|I becomes 

is ps k|| — (asin0cos#)k^ — ► kii, (B3) 

if ps (a sin 9 cos 0)k|| + k ± -> k ± . (B4) 
In this limit, ^ s is mostly proportional to k|| and ^/ to kj_. When a — » oo, equations l|B2(l and (|B 1 f) becomes 

^ s « -0+ (sin6»cos6»/a)k, (B5) 

|/ PS (sin<9cos6>/a)(9 + k. (B6) 
When a = oo, slow modes are called pseudo- Alfvenic modes. 

The quantity V • v; becomes ~ (vi/l)(k ■ ^), where (...) denotes average over Fourier components that have k ~ l/l. 
For slow modes in low-/? plasmas, k ■ £ s ~ k • k|| ~ kn/k and, therefore, 

V • Vl^low ~ Ui.sioui/Zn. (B7) 

For slow modes in high-/? plasmas, k • £ s ~ k • (sin6'cos6'/Q!)k ~ (ku/k)/a and, therefore, 

V ' ^l.slow ~ v Lslow /(/3Z||). (B8) 

For fast modes 

V -Vijast ~ Vljast/l (B9) 

for all values of /?. Since Zii > Z, the value in eq. (|B9|I is larger than that in l|B7() . This explains why fast modes are 
more efficient in particle acceleration. 

Estimates for density fluctuations 

In CL03, we provided estimates for density fluctuations in MHD turbulence. In case there are only slow modes, the 
density fluctuations are 

(Sp/p) s i ow ~ VL.siow/a (BIO) 

in low-/? plasmas (CL03). However, when Alfven modes are pres ent, slow modes are enslaved by Alfven modes (see 
discussions in Lithwick & Goldreich 2001). In this case, eq. (|B10|) is modified to 

(Sp/p) ~ V L , s low/vA, (Bll) 

which is ~ 0(1) when vl >s i ow ~ va- The contribution from fast modes is {Sp/p)f as t ~ VLJast/vA- Therefore, when 
VLjast < vl.sI ow, slow modes dominate density fluctuations in low-/3 plasmas. Although fast modes can compress 
more (c.f. eq. (|B9|) L short time scale, or wave period, of fast modes makes them not so efficient. 
In high-/? plasmas, fast modes dominate density fluctuations: 

(Sp/p) ~ v LJast /a. (B12) 

When VLjast ~ va, we have (Sp/p) ~ l/yf]3. 
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